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The temperature inversion properties of the internal energy, E, on odd spheres, 
and its derivatives, together with their expression in elliptic terms, as expounded 
in previous papers, are extended to the integrals of E, thence making contact with 
the theory of modular forms with rational period functions. 

I point out that the period functions of (holomorphic) Eisenstein series com- 
puted by Zagier were already available since the time of Ramanujan and I give a 
rederivation by contour integration. Removing both the Planck and Casimir terms 
gives a fully subtracted form of the series which allows a more elegant and compact 
treatment. I expound the relation to Eichler co homo logy cocycles and also rewrite 
the theory in a distributional, Green function way. 

Some historical and technical developments of the Selberg-Chowla formula are 
presented, and it is suggested that this be renamed the Epstein-Kober formula. 
On another point of historical justice, the work of Koshliakov on Dirichlet series is 
reprised. A representation of a 'massive' generalised Dirchlet series due to Berndt 
is also reproved, applied to the Epstein series and to a derivation of the standard 
statistical mode sum, interpreted as a Kronecker limit formula. 
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1. Introduction. 



In two previous works, [1], [2], we have investigated the thermal quantities for 
scalar and spinor fields on the space-time TxS d , where d is odd. The conformal 
scalar internal energy, E, for example, is a linear combination of 'partial' energies, 
e t , t = 1,2,..., (see below). These quantities possess a known behaviour under 
modular transformations and are central to elliptic function theory. In particular, 
a temperature inversion symmetry can be exhibited. 

From basic elliptic properties, it was shown that E could be expressed as a 
polynomial in e 2 and e 3 , while, for the specific heat, and all higher derivatives, one 
has to include e 1 , preferably via a modular covariant derivative. The free energy 
involves an integration and things are not so simple as obstructions arise to simple 
modular behaviour. It is this aspect that is explored in the present work. I use this 
thermal angle just to motivate the introduction of various quantities, the analysis 
of which allows us to make contact with various topics in the theory of modular 
forms. Some of the procedures can be given a physical terminology, which might 
be suggestive. 



2. Mellin transforms and the period polynomial. 

It is advantageous to begin from the basic Mellin transform, used e.g. by 
Malurkar and Hardy, for the Eisenstein series (or partial internal energy, [2] [1]), 

B 2t , ~ n 2 '-y™ 



q 2n 

»=i (1 ) 

j j rc+ico 

= -( R (l-2t) + — J dsT{s)Q R {s)Q R {s-2t + l)(2nb)- s 1 

where q = e l7rr and c lies above 2t. I am using b = —ir = (3/2w = l/£ as yet 
another convenient variable. Although, physically, b is real, if we wish to consider 
general modular transformations then we must allow it to become complex. 

The Mellin transform has been used systematically in finite temperature field 
theory and elsewhere in discussions of asymptotic limits. The survey by Elizalde et 
al, [3], is useful and I further mention, as being somewhat relevant here, Cardy [4], 
Kutasov and Larsen, [5]. Terras, [6] pp. 55, 229, can be consulted for aspects of the 
mathematical side. In number theory see Hardy and Littlewood, [7], and Landau, 
[8,9]. 
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Equation (1) can be looked upon as a continuous expansion in b. Various 
limits in b are uncovered by displacing the integration contour. The pole in the 
integrand at s = 2t corresponds to the Planck term, while the remaining one at 
s = corresponds to minus the Casimir value, which is the first term in (1). The 
high temperature limit (b — > 0) follows by pushing the contour to the left, beyond 
the Casimir pole, the contribution from which cancels the first term in (1) leaving 
just the effect of the Planck pole and a correction that tends to zero exponentially. 
Conversely, on moving the contour all the way to the right only the zero temperature 
Casimir term remains. 

This cosmetic asymmetry can be avoided (if desired) by also extracting the 
Planck term, defining, 

er b (b)=e t (b)-ic«(l-2t)(l + (zb)- 2 *), (2) 

and writing 

e " b ( b ) = ^7 / dsT(s) ( R (s) ( R (s -2t + l) (2nb)-\ (3) 

where now c, at the moment, lies anywhere between and 2t. 

The quantity in the integrand is the Hecke L-series of the normalised Eisenstein 
modular form, 

L(G t ,a) = C R (a)C R (a-2t+l), (4) 
satisfying the (typical) functional relation 

(27r)- s r( S ) L(G t , s) = (-lY(2ny- 2t T(2t - s) L(G t , 2t - s) . (5) 

I now consider a multiple integration with respect to b of the internal energy. 
Actually for the free energy only one integration is needed but the general case is 
mathematically important. So integrate e^ ub , h times to give, using (3), 

m Jx (6) 

with h < c < 2t. 

I have imposed the boundary condition that each successive integral vanishes 
at x = oo, i.e. at zero temperature. 
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More poles have been introduced into the integrand from the F function and 
the contour has to be moved to the right to avoid them but we should stop when 
they pass beyond the contour at Re s = c = 2t — 5, < 5 < 1. Hence the maximum 
h can be is 2t — 1. 

I now set h equal to this maximum, 2t — 1, simply for the reason that it gives 
the ultimately recognisable, and important, quantity, 

\2t 



& t (x) - ^|^T) / x °° db ( b - x ) 2t " 2 e r>) 



2 r 



C+ioo p/ _ rw , i \ (^) 

■»» ( 2j ),-2, ' <«(') C,(» - 2t + 1) x-+«-' 



with 2t — 1 < c < 2t. The first line is recognised as a Weyl fractional integral. 

The significance of this quantity for statistical mechanics on spheres is admit- 
tedly obscure because the internal energy is a sum of the e t for different t, and the 
integral depends on t. Only for the 3-sphere, where t = 2 and E = e 2 , is there a di- 
rect relation but even then 4 is a third integral of E having no particular thermal 
meaning. Nevertheless, its further analysis is not without interest as it connects 
with some salient mathematical concepts. 

According to my general programme, interest lies in the behaviour under mod- 
ular transformations. My treatment is equivalent to that of Malurkar, [10], c/also 
Guinand [11]. In contrast to the derivatives of E, there are obstructions to the 
modular invariance of the integrals and the period polynomials or, rather, period 
functions, provide a measure of these obstructions, as we will see. 

Firstly under translations, b — > b — i, and, because of the Planck subtraction, 

er b (b - i) - er b (b) = (-Y +ll -( R (l - 2t) {j^r t - -L) , (8) 
which implies, through (7), that 

2t (x-z)-0 2t (x) = 2z-^£^. (9) 

XIX t J 

Under inversion, x — > 1/x, the integrand involves 

er b (l/b) = (-l)V*er>) 
since the regularising subtraction, (2), maintains the inversion property, 

e,(l/b) = (-l)VS(b), 
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enjoyed by the full Eisenstein series in (1), 
Therefore, 



Ml/x) =^1^ f x (b - 1/x) 2 *- 2 er b (b) 



2(2tt) 



L/X 
2t /-0 j k / 



/ x ^(l/b'-l/x) 24 - 2 ^^!/^) (10) 



T(2t-l)J x 

=(-1)' x 2 " 2t T^iy jf rfb ( x - b ) 2t " 2 e r b (b) 



The amount by which 2t violates the pure modular inversion property is the 
quantity, 

^2t(x)-^x) 2t - 2 ? 2t (l/x)=P,(x), (11) 

or, from (7) and (10) 

^ (x) =r§r-T) [ db (x ~ b)2 " 2 erb(b) ' (12) 

The calculation of P t from the Mellin contour form goes as follows. From (7), 
&t(l/x) 



,2t-2 



2m J 2 t-c-ioo (2tt) 



by setting s — > 2t — s. 

The ^-functional equation (one can also, more easily, use (5)) leads to, 

( R (s)( R (s-2t + l) 

= 2 ply-° sin ( ws W C«(i - «) 2 ( r 2 ( 7r 2 ) t 2 7_! ) - 2i + i)/2)C B (2« - «) 

- " 1 ' ' S)l * 2i " " sin(7rs/2) sin( ;;(>• - 2t + l)/2)C fl (l - s)Ci?(2t - s) 



(2n) 2t + 1 ~ 2s 

= (-i) t 2r( [ 2 ;) a 2+i? a r a) - -)c h (m - -) 

(-l^rci-s 



(27r) 2t - 2s r(s - 2* + 1) 



C fl (l - s)C fl (2t - s) , 



and so (13) becomes, 

02t(V><) 

2 (- 1 ) t " 1 1 f 2t - c+t °°, r(a-2f + lK . ^ (14) 



X 2 *" 2 27TZ 



J2t-c-ioo \ zn ) 



I then find the alternative expression for the obstruction, P t (x), defined in (11) or 
(12), 



pm-U r°°- [ °°\dc T{s - 2t + 1) c (s)c r a -2t+nx»— 1 

Pt(X) "^Uc-*oo L-c-ioo ) (27T)-« Cfi(S) ^ (S 2t + 1)X 

Translating the contours so that they mutually cancel leaves the contributions of 
the intervening poles i.e. those between and 2t. Thus 

= i / r( ?~?_t 1} C*00 Cn(s -2t + l) x 2 — 1 



7T* J c (27t) s " 2 * 

= (-1)*" 1 - / g^Cr^Cr^ - s) cosec (tts/2)x 2 *- s - 1 
i Jc 

t-i 

= 2n( R (2t - l)((zx) 2 *" 2 - 1) - Mj2C R m C R (2t - 2j)(ixf t - 2 ^ 

3 = 1 

= 2n( R (2t - !)((*,)»-» - 1) - K-W ^^ (fa)*-"-' , 

(15) 

which is a polynomial of degree (2t — 2) and is the final, main result of this section. 
The odd powers come from the poles of the cosec while the two even powers arise 
from the poles of the Riemann ^-functions. P t (x) is real if x is. 

Setting x to zero in (15) and (12) produces the specific formulae, 

2(27T) 2t 1 f°° ,2t-2 sub, 



( R (2t-1) = -/ dbb 2t - 2 tT b (b) 

R } r(2f-l) J 1 1 } 

using inversion. These can be checked numerically. More generally, expanding (12) 
gives the other nonzero moments of e| ub (b), 



dbb^-\r h (b) = ^^- ) B 2j B 2t _ 2j , Kj<t-1. (17) 
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The periods given by Kohnen and Zagier, [12], are identical to (16) and (17) and I 
now have made contact with this mathematical notion. 

The periods of a cusp modular form, JF, of weight 2t are the numbers (e.g. 
Lang [13]), 2 

/■OO 

r n (T)= dbF{ib)b n , < n < 2t - 2 . (18) 
Jo 

The end points are cusps where T vanishes (exponentially) to ensure convergence. 
The restriction on n is a polynomial one (see [13].) 

It is important for the following to assemble the periods into a period polyno- 
mial, which, as usual, is a handier quantity, 



n=0 ^ ' 



(19) 

/■OO 

= / db(x- b) 2t ~ 2 F(ib). 
Jo 

In order to extend the notion to non-cusp forms (such as Eisenstein series) one 
course (see Zagier, [14] and Kohnen and Zagier, [12] §4), is to continue the integer, 
n, into the complex plane by employing the Hecke L-series of JF via the Mellin 
transform, 

L{T, s) = ^ r dt (F(it) - F(ioo)) t*- 1 , (20) 
1 K 8 ) Jo 

which converges for s > 2t — 2 and, continued in s, satisfies the functional relation, 

(2tt)- s r(s) L(F, s) = (-l) t (2n) s - 2t F(2t - s) L(F, 2t - s) , (21) 

by virtue of the (inversion) modularity of J 7 , and conversely. This is usually at- 
tributed to Hecke, [15] (but see the Appendix). Terras, [6] p. 229, gives a useful 
pedagogical treatment of this topic. 

Then the periods could be defined, in general, by, [12], [14], 

r ^= lim + i ^S L ^ , ^ 

s^n+i (2tt) s ^22) 

= lim L*(jF, s) , 

s— >n+l 

because this coincides with the integral expression, (18), for cusp forms. 



I prefer the normalisation of Kohnen and Zagier, [12]. 
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The explicit values of r n (n = 0, . . . 2t — 2) for Eisenstein series (which is all 
that is required because of a decomposition theorem for modular form space) are 
easily computed from (4) and are given in [12]. They agree with (16) and (17), as 
already noted. 

However, the Dirichlet series, L*(jF, s) has poles at s = and s = 2t. A means 
of incorporating these values is given by Zagier, [14], using the period polynomial. 
One notes that the summation in the cusp period polynomials can be extended, 



nE Z 

= E(- 1 ) n ( 2 ^ 2 )^ n+lL *(^^+ 1 ) x: 

n€ Z 



(23) 



because the binomial coefficient vanishes outside the range n = 0, . . . , 2t — 2. For 
non-cusp forms, however, the sum extends from — 1 to n = 2t — 1, the nonzero end 
values arising from the singularities of L*{T,s) mentioned above. The resulting 
expression can then be taken as a definition of the period polynomial of a non-cusp 
form. 

Again, for the Eisenstein series, G 2t , = tf, the computation is straightforward, 
[14], and (23) with (4) yields, (n = 2j - 1), 

Tg{x) = f^F ( 27rC « (2t " 1} ((" 1 ) t " lx2t " 2 " !)" 

B B \ ( 24 ) 

(27F) ^ !) m {2t-W )' 

which is 1/x times a polynomial. For comparison with [14], Zagier's X equals my 
zx. 

I now comment on this construction. The periods, and period polynomials, for 
non-cusp forms are defined by analogy to those of cusp forms, the justification being 
that the expressions, for jF(ioo) equal to zero, reduce to those for cusp forms. The 
means to attain this is not unique. Instead of the 'zeta-function' regularisation, as 
in the definition of L, (20), I employ the subtracted form, as in (2), 

^ sub = ^-a -^, a = ^(zoo), 

whose moments are non-infinite and provide an alternative definition of the periods 
of a non-cusp form which equally reduces to that for cusp forms when a = 0. This 
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subtraction maintains inversion modularity, but violates translational, whereas the 
subtraction of just a does the reverse. Making use of my earlier considerations of 
Eisenstein series, see (12), I define the subtracted period polynomial as an integral, 

poo 

r^ b (x) = / rfb(b-x) 2 *- 2 ^ sub (b). 
Jo 

For Eisenstein series, (12) shows that r£ b (x) is proportional to P t (x) whose 
explicit form is given in (15). 

It is necessary to compare this definition with that of Zagier, (24). The dif- 
ference is the absence of the 'end point' terms proportional to 1/x and x 2t_1 . This 
can be remedied by including, in the terminology of this paper, the regularised con- 
tributions from the subtracted terms, the Casimir and Planck terms, which can be 
accomplished simply by expanding the closed contour in (15) to include the corre- 
sponding poles at s = and s = 2t. This just extends the sum by the two end 
points giving the well known quantity, 

t 

R t (x) = 2n ( R (2t - 1) ((zx) 2 *" 2 - l) - 4i £ C R (2j) ( R (2t - 2j) (zx) 2 ^" 1 , (25) 

3=0 

which agrees, up to a factor, with Zagier's expression given in (24). 

Furthermore, the use of the subtracted form, jF sub , allows one to give a cleaner 
integral expression for the enlarged period polynomial than that in [14], i.e. , 

r -i(x) = J™ db (b - x) 2t " 2 ^ sub (b) + (x 2 *" 1 + x" 1 ) . 

While equivalent to the method in [14], I believe the contour approach is neater 
and is capable of being extended to the general theory of periods. Zagier, [14], shows 
that logistic simplifications occur in this theory when it is extended to include non- 
cusp forms such as the Eisenstein series. 

As an historical point, it will be noticed that the above expressions have been 
available since the time of Ramanujan, with contour derivations by Malurkar, [10], 
and Guinand, [11]. In order to investigate this point further, I return to the mode 
sum expressions for the thermodynamic related quantities. 
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3. Back to summations. 



I start with the Eisenstein series forms and work the development to parallel 
the quantities of the previous section. This allows me to introduce the useful cal- 
culations of Smart, [16], on the Epstein function (which will arise later) using its 
relation with Eisenstein series. 

The subtracted partial energy ef ub in (2) corresponds to the the 'regularised' 
Eisenstein series, 

oo oo 

G r h (r) = y' y' - — 

t v ; mr + n 2 * (26) 

m= — oo n= — oo v ' v ' 

= G t (T)-2( R (2t)(l + T- 2i ), 

The relation is, [2], 

t sub/ \ (2t - 1)! —sub 

€ < (T)= 2(2^ - Gt (r) ' (27) 
Now, corresponding to (7), integrate (26) with respect to r from r to zoo, 2t — 1 
times to give, formally, 

f^TyMr), (28) 
where 4> is defined by, [16], eqn.(2.11), 

oo oo 

y — o7 — r • (29) 

m= — oo n= — oo 

The behaviour of the series 4>2k( T ) un der modular transformations, in particular 
r — > — 1/t, has been especially considered by Smart, [16], who refers to the 0's as 
'modular forms of weight 2 — 2t with rational period functions'. 

Using the partial fraction identity (see also Hurwitz, [17], Eisenstein, [18]), 

1 1 r 2 *" 1 1 ^ T 2t-J-1 

m zz i mr _|_ n n zi l mr _|_ n m j n zi j 

J = l 

one finds easily the inversion relation, [16], eqn.( 1.10b), 

t-i 

cj> 2t (r) - r 2t - 2 4> 2t {-l/r) = -4 £ r 2 *- j y- 1 C fl (2j) C fl (2t - 2j) . (31) 



Also, under translations, 



2t (r+l)-0 2t (r) = 2-^£^. (32) 
10 



Up to a simple factor, one might expect 4> 2t to be the same as 4> 2t in (7). 
However the cocycle functions in (31) and (15) differ by the two even powers. I 
note that these are related to the j = 1 and j = 2t — 1 terms in the sum in 
(30) which have gone out in the passage to the difference, (31), by antisymmetry. 
However the summation over m, respectively n, for these values of j is conditionally 
convergent and is therefore subject to ambiguity, that is to say, a choice has to be 
made, as in the discussion of the Eisenstein series, G 2 . The regularised Mellin 
transform approach results in a different definition of the series (29) to that used by 
Smart, [16], p. 3. The relation is easily found by taking the even powers in (15) over 
to the left and then we see that one has the relation, (further remarks are given 
later), 

02t(x) = i0 2t (r) - 2nQ R (2t - 1) , x = -ir , (33) 

also agreeing with (9) and (32). 

Smart, [16], also gives the relation with a known Lambert series as 



SJt) = Y r - „ = — — <Mr) — ( R (2t) - -( R (2t - 1) 

(34) 



m=l 



4tt vrziv ' x 
= 4^^ (X) ' 

where is the function having i? 2i (x), (25), as its inversion cocycle function, 

^2t(x) - (zx) 2 *- 2 ^ 2t (l/x) = R 2t (x) , (35) 

as is easily confirmed. 

The introduction of ip has restored translational invariance, (the Planck term 
has been put back in (34)), 

^ 2t (x-i)-^2tW = 0. (36) 

It is almost convention to consider modular integrals, for example, to be periodic 
(e.g. Knopp, [19]). See some later comments. 

Equations (35) and (34) imply an inversion relation for the series S t which has 
a certain history, some of which is detailed by Berndt, [20], p. 153, who says that 
it was first written down by Ramanujan. As is clear, it is quite the same as the 
expression for the Eisenstein period polynomials. 
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The Lambert series, S t (r), has been considered by Apostol, [21], Berndt, [20], 
Grosswald, [22], and Smart [16], for example. It can be related to statistical me- 
chanics in the following way. 

One can resum the statistical expression for the partial free energy in a standard 
manner, beginning, (q = e KlT = e _7r /^), 

j. oo £ oo 2t — 2 

ft = - i E " 2 - 2 mi - f 2 "> = ««. - 1 E V' 2 "" ■ (37) 

n=l m,n=l 

The scaled entropy is defined by, 

\ ^> / m,n=l 

Concentrate now on the summation and rewrite it, 

E V" 2 "'" = g2 " 2 E ■ (39) 



m,n=l 



— ' m 2t 1 1 — q 2 

1=1 



in terms of S t (r), (34). 

This Lambert series has been discussed in connection with generalised Dedekind 
ry-functions by Berndt, [20], and it is worth noting that, in a certain sense, log?7 is a 
modular form of zero weight, in agreement with the results e 1 = —D log rj, [1], which 
could therefore be written, e 1 = — Vlogrj, in terms of the covariant derivative, V. 

The question, not answered here, is whether there is an 'expansion' or repre- 
sentation theorem for modular integrals combining that for modular forms (of any 
real weight, see [23]). It is not likely that one can find a pure elliptic formula. As 
evidence I cite the particular evaluations of Smart who gives at the lemniscate point 
( T = i,x=l), 

? 4 (1) = ^-27^(3), 
and, more generally, if 2t = 0(mod4), 

2tt 2 ' 2tt 



s t (i) = ^ - Cj ^f^ + ^ B-i) J+1 c*(2t - v)t R m 



a result actually due to Lerch which follows from (35) with (25). 

As a final general comment in this section, it is always possible to obtain the 
required quantities as (infinite) (/-series by direct integration, but these do not count 
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as closed forms. For example, one integration yields the (/-series for the partial free 
energy in terms of the arithmetic functions, a k (m), (the sum of the k-th powers of 
the divisors of m), 

I 00 q 2m 

m=l 

which amounts only to doing the n-summation in (37), or an integration of the 
standard relation, (1), 

e t = -^+£*2t-iM<Z 2m . 

m=l 

One also has the formula, e.g. [24,25], 

m=l m=l 

in agreement with the result of a (2t — l)-fold integration. 



4. Another approach via modular integrals. 



The period functions are cocycle functions associated with Eichler cohomol- 
ogy and modular integrals. These notions put a different slant on the preceding 
formulae. In effect one begins again. 

The best place to start is the important (but particular) analytical result of 
Bol, [26], which states that, for any reasonable function, <p(r), 

(D r+1 <p) ( 7 r) = (cr + d) r+2 D r+1 ((cr + df ^( 7 t)) , (41) 

where r is a nonnegative number, and I use, 

(^ + V)(r)=(^) r+ V(r), 

with, 

7 =(* c *)eSL(2,R). 

This shows that, if ip is an automorphic form of weight — r, then D r+1 tp is an 
automorphic form of weight 2 + r, see also Petersson, [27]. Note that it is the 
ordinary derivative that enters here. The proof of (41) follows by induction. 
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The classic discussion of modular integrals can proceed by asking for the inverse 
of BoPs result. That is, given a form of weight 2 + r, p(r), as a source, can one 
integrate the differential equation, 

D r+ V = P (42) 

to obtain (p, a form of weight — r? Clearly this will not be so in general because of 
the freedom introduced by constants of integration. That is, one would expect <p to 
behave as a form of weight — r up to an additional piece, the general form of which 
can be determined from the differential equation (42) and Bol's theorem, (41). 
Generally we would write (ignoring any multiplier systems for simplicity), 

rtyr) = (or + d)- r (<p(r) + P( 7 ,r)) , (43) 

where -P( 7 , r) has to be a polynomial in r of order < r, The inverse to Bol's result 
is that, if any <p satisfies (43), then p, of (42), is a form of weight (2 + r). This 
follows most simply by substituting (43) into (41). 

P must satisfy the 1-cocycle consistency condition, [28], 

P{lil2, r) = (c 2 r + d 2 ) r P( 7l , 72 r) + P( 72 , r) 
^(7i,r)| 72 + P(72,r). 

(p is known as an automorphic or Eichler integral, of weight — r, and P is the 
associated period polynomial, [28,29]. I have also introduced the useful, standard 
'stroke' operator, | 7 , sometimes denoted |[ 7 ]. 

The general solution (equivalently the indefinite integral) of the differential 
equation (42) is a particular integral plus the complementary function (a zero mode) 

V {t) = ^-^y I da (r - a) r p(a) + 0(r) , (45) 

where 0(r) is a polynomial of degree r. The integral, in the upper half plane, is 
path independent if p has weakish analytic properties and, if desired, the solution 
can be made independent of the lower limit by suitable adjustment of 0(r). In 
elementary fashion, 0(r) is determined by the first r derivatives of (p(r) at r . 

All this is general analysis. In modular applications to SL(2,Z), r is chosen 
as the cusp ioo and the limits in (45) reversed. In terms of Fourier expansions, a 
holomorphic cusp form source is 

P(r)= $> m ? 2m , ( 4 6) 

m=l 
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and the particular Eichler integral of p, is, from (45), 

oo 

^) = £^ 2m > ( 47 ) 

m=l 

which has the same periodicity as p, a result of choosing r = zoo. 

The problem in the present Eisenstein case has already been noted. The Eisen- 
stein series, G t (r) does not vanish at zoo (or 0). It is a 'non-cusp' form and the 
direct integration of (42) meets obstacles at an elementary level. This can be appre- 
ciated simply from the appearance of 1/m factors during integration, and m can be 
zero in G t . This fits in because the m = terms in G t are those making G t nonzero 
at too, the zero temperature value. Saying it yet again, the Fourier transform of G t 
begins with an m = term, the Casimir contribution, so (47) makes no immediate 
sense. 

One way of avoiding these obstructions, is to regularise G t by dropping the 
offending term(s) as in (26) and consider the differential equation, 

D 2 '- 1 <P 2t = -4^ ub , (48) 
which has the particular integral, <p 2 t = 4>2t where, 

(9 -\2t-l rioo 

0»(r) = / dr' - r) 2t ~ 2 er b (r') , 

the constants and boundary conditions being chosen so as to reproduce (7). The 
upper limit is a cusp. Smart's function, <f> 2t , results if one adds a complementary 
constant (and multiplies by i). (See (33).) 

Because of the subtraction, the quantity ej ub (r) is not exactly a modular form, 
rather one has, expressed slightly generally, (r = 2t — 2), 

e ( 7 r) = (cr + d) r+2 (e(r) + E(>y, r)) (49) 

and also 

0( 7T ) = ( CT + d)- r (<P(r) + P( 7j r)) , (50) 
which are related by the same differential equation, (42), 

D r+1 (j)= -47T€. 

Bol's formula, (41), relates P and P, 

D r+1 P = -4nE, (51) 
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and these relations generalise the usual ones applicable to cusp forms. The co cycle 
formula (44) still holds of course, as (50) is the same as (43). 

I remark that <j> and e sub , not being periodic, do not quite possess Fourier 
expansions, as in (46) and (47). However ip, in view of (36), does, 

oo 

?(x) = E^^ 2m ' 

m=l 

Because P, this time, satisfies (51), it is not required to be a polynomial. 
Generally, since the modular group is generated 3 by S and T, P(7, r) is a rational 
period function. In the present case, however, P(S, r) is a polynomial, given by 
(15), 

P(S,r) = P t (r), (52) 
while the (non-zero) P(T,r), which follows from (9), 

P(T,r) = 2C R (2t)(^--^- j y (53) 

is a ratio of polynomials. Of course, P(±l, r) = and I remark that the polynomial 
P(S,t) runs from r° to r 2t ~ 2 . 

As a simple check of the algebra, set r ) 1 = 7 2 = S in the cocycle relation, (44). 
This gives, 

P(S, t)\(1 + S) = P(S, t) + r r P(S, St) = P(S, t) + r r P(S, -1/r) = 

which, for r = 2t — 2, agrees with the explicit form, (15). 
Also, as examples, 

P(TS,t) =2( R (2t)— + P t (T) 
P(ST,T)=l^ )+ P t (l + T). 

I note that the period function appears to have poles at —1, 0, +1 and zoo. 
From (49) and the form of e sub , (2), one finds 

E(S,r) = 0, E(T,t)=( j -^--^) 

3 Beware that some authors, e.g. Knopp, swap the usage of the symbols S and T. 
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which are consistent with (53) and (52). 

For completeness, and to aid comparison, this is a convenient place to outline 
some standard elementary things about period polynomials, mostly for SL(2,Z). 
Two basic 'polynomials' are the inversion one, P(S,r), and the translation one, 
P(T,t). For cusp forms, the latter vanishes because of periodicity, P(T,t) = 0, 
[19], corresponding to the periodicity of the modular integrals, c/(46), (47). In this 
case, P(S, t) satisfies the Eichler-Shimura relations 



with 5" 2 = 1 and (TS) S = 1. A standard way of proving these uses the complex 
integral form for P, e.g. [14], but they can also follow, algebraically from the cocycle 
relation as I now show. (See also Knopp, [30], for slightly different details.) 

For the first identity, just set, as above, 7i = 1 and 7 2 = S in (44) and it falls 
straight out. The second identity is slightly more work. Using (44) twice, one has, 
dropping the r argument temporarily and noting (TS) 2 = (TS 1 ) -1 , 



P(S)\(1+TS + (TS) 2 ) = 

P(S) + P(STS) - P(TS) + PiSS^T' 1 ) - P(TSTS) . 
Further, setting r ) 1 = T in the cocycle condition, and using the assumption 
that P(T) = 0, yields 



One then sees that the first and third terms on the right-hand side of (55) cancel, 
as do the second and fifth. The fourth term vanishes on its own because (56) shows 
that P(T _1 ) = 0, and (54) follows. 

In one development, the notion of modular integral has been somewhat di- 
vorced from that of integration. One says that a modular integral is any function 
(with appropriate analyticity properties) that obeys the quasi-modular relation 
(43), where the period function P must still, of course, satisfy the cocycle condi- 
tion (44). Motivation is then provided to look for modular integrals with rational 
period functions, i.e. ratios of polynomials. In this case it seems to be conven- 
tional to assume that the translational period function, P(T,t), vanishes, [19], so 
that both Eichler-Shimura relations, (54), still hold. There is a certain interest in 
computing such rational period functions and investigating their zeros and poles. 
The easiest introduction to this topic is by Knopp, [30]. It is, however, somewhat 
peripheral to our considerations which are more concerned with non-cusp forms i. e. 
the Eisenstein series, e.g. [14]. 



P(S,t)\(1 + S) = 
P{S, t)\(1 + TS + (TS) 2 )=0 



(54) 



P(T 7 ) = P( 7 ) . 



(56) 
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5. Green function distributional description. 



Developing the formalism a little further, it is possible to give the discussion 
of Eichler integrals a distributional or Green function look. For this purpose, it is 
better to use the 'real' variable form and consider the differential equation, 

£^^(x) = p(x), (57) 
for x > with solution vanishing at oo, cf (7), 

/•OO 

#c)= / dx'$ 2t _ 1 (x'-x)p(x') (58) 
No- 
where the generalised function, $ Q (x) , is 

*«(>) = £• (59) 

The generalised function is concentrated on the positive x-axis, i. e. is equal 
to x a for x > and is zero for x < 0, e.g. Gelfand and Shilov, [31] I, §5. 5. $ 2 t-i 
acts as a Green function satisfying, 

^ T $ 2t _ 1 (x'-x)=5(x'-x). 

The formal, distributional equivalent is obtained from the convolution, valid for all 

a and /?, 

$a * $/3 = $ a +/3 , (60) 

by setting a = —(3 = 2t — 1 and noting 

*_ fc (x) = <5<*>(x), k = 0,1,2,.... (61) 
An important formula is the Fourier transform of <E>, [31] I. p. 360, 

and the inverse, which is a (continuous) eigenfunction expansion, 

*(x) = j™/ 2 ±- / rfa-^ -. (63) 
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In the present case, the source p is assumed to be periodic under x — > x — i (the 
Planck term is not removed) and to converge as x — > oo. So one has the 'Fourier' 
series, 

oo 

P(x) = $>- e " 2m?rX ' x>0 - ( 64 ) 

m=l 

Substitution of this and (63) into (58) gives, 

^(x) = (-l) 1+1 -y^ da- -. (65) 

Closing the contour in the upper half a-plane I obtain, x > 0, 

1 oo 

*m=mhE^^ x - («> 

which, it is no surprise to see, is just (47) derived by a much longer route but one 
on which we encounter some handy information. It is seen also that the source 
periodicity has been transmitted to the solution. 

Including an m = term in the Fourier series, (64), still gives convergence at in- 
finity, but of course the solution, (66), breaks down, as has already been mentioned. 
As with all Green function approaches, when this happens the corresponding eigen- 
function has to be removed from the source. This can be done formally by redefining 
the Green function to exclude this function which is, in the present constant 
zero mode. Altering the contour in the Fourier integral (63) gives a modified Green 
function, 

-i i>oo+i8 p —iaX 

*"£ l( x) = (-!)' — / rfa^, 0<5<2, (67) 

Z7U J-oo+i5 a 

which automatically takes out any constant part of the source. This can be seen by 
going back to (65) and replacing the integral by one from — oo + id to oo + iS, 

1 °° i-oo+iS iaX 

m = _ £ ft „ *, a2t _ 1(a _ 2mn) . (68) 

For < 5 < 2, if x > 0, closing the contour in the upper half plane still yields (66). 
However now, even if the source Fourier series, (64), is extended to the constant 
term, m = 0, the result is again (66). This should be denoted by ip suh (x) and one 
can write _ 

/"OO 

^ sub (x)= / dx'^ i (x'-x)p(x'). (69) 
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Retaining the general structure, (63), for any cm, one can define for any a, 



^ roo+iS g—iaX 



$ sub (x) = e «ra/2 / da——. (70) 



2tt J-vo+is a c 

Then $^ ub (x) is still concentrated on the positive x-axis and the convolution (60) 
also remains valid, 

K uh * <^ ub = K% (7i) 

together with (61), 

To check these statements, firstly, if x < 0, closing the contour in (70) in the 
upper half plane yields zero as the only singularity is at a = 0. Then 



W b „ $ s ub) (x) = / ^ / d l da> 



where the lower limit on the x' integral is allowed because $^ ub (x ; ) is concentrated 
on the positive x' axis. 

Setting a = s + id, a' = s' + id, I get, in detail, 

i poo ,>oo ,>oo -iaX-iX' (s' -s) 

-iaX 



_ e iir(a+f3)/2 



_ e i7r(a+/?)/2 



^ i>oo poo g — wtX 

— / ds ds' t^S(s~ s') 

1 r°° p -iaX i roo+io 

±- / ds^=e^ a +^ 2 ± da- A , 



Equation (72) follows either from (71) or (70) since the integrand has no sin- 
gularity, and so one can set 5 = 0. 

Eichler cohomology. 

As a side remark, the general solution to (57) consists of the particular integral, 
(58), plus a solution of the homogeneous equation, i.e. a zero mode, which in this 
case is a polynomial, 0(x), of degree at most 2t — 2. The contribution of this to 
the cocycle polynomial, P of (43) , constitutes an element of the set of coboundaries 
and a cohomology of polynomials can be set up (Eichler cohomology). Working 
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with the cohomology classes removes the constants of integration ambiguity in the 
solution of (57). That is, there is a one-to-one correspondence between p and an 
element of the cohomology group denoted i? 1 (r, II r ), where Ii r is the vector space 
of polynomials of degree < r = 2t — 2 (for cusp forms) and V is (here) the modular 
group. The development of these ideas, while significant generally, is not required 
here. 



6. Epstein approach to thermal quantities. 

An 'earlier' form of the free energy, or effective action, follows from the thermal 
^-function, which, in this case is related to the Epstein ^-function. Starting from 
this I will be able to make contact with the Eisenstein formulation and derive 
useful relations. Using the Epstein function is equivalent to the non-holomorphic 
Eisenstein series. As a rule the analysis is harder. In reality one is doing twice the 
necessary work. 

In some ways we are going backwards. The use of holomorphic forms, e.g. (27), 
should be considered a simplification available through working with an explicit 
square root of the Laplacian. 

Again splitting up the degeneracy, one is led to define a 'partial' ^-function, 

OO 9+_9 

Ct(S,/3) = 20 £ (4ttW//? 2 +n 2 /a 2 )* ' ( ?3 ) 

— OO 

and the related free energy, 



oo 



F 



-Lu m lY'-J?— — . (74) 



2t-2 



t 



I have rescaled the free energy by defining, as above and as in [2], F t = aF t 
and have also taken a preliminary limit of s — > in an extracted factor of (2tt/(3) 2s . 
This will be justified in a moment. 



F-( d V" 1 lim r(s ~ t + 1) V' 1 f75l 

Not surprisingly, the degeneracy factor is replaced by a derivative and the 
resulting sum is precisely an Epstein ^-function. This has been used before e.g. 
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[32]. The point is that there is no pole at s = 0, which justifies the limit. In terms 
of the Epstein function, Z 2 , 



lW 5 / d r(«-t + i) 



F ' = ( - 1) ^UvaJ il m „%nr Z2(2s - 2i+2 ' A) ' (76) 

where A is the matrix (the 'modulus'), 

1 

o r J 

The idea now is to use the standard functional relation satisfied by Z 2 (s,A), re- 
peated here for convenience, 

2s-l p/i _ \ 

Z 2 (2s, A) = — f K 4r, ^ 2 (2 - 2a, A~ l ) , (77) 

to give 

P _L ( _ —4—1 * ^ lim n 2s ~ 2t+1 F(t ~ ^ Z 2 (2t - 2a, A" 1 ) . 
1 8tv\ d(l/£ 2 )) S->o r(a + l) 2V ' ; 

Z 2 (2s, A) has a pole at a = 1 only, so working at t > 1 to begin with, one has 



F -t tW^EW f d \~\ Z (2f A~ x \ 



(78) 



As the simplest case, set t = 2 giving the three-sphere. Then 



F = ^ ^Wm^ A ^ = "i^l f (79) 

agreeing with the result in [32] , allowing for the changes in notation. 

This formula can be taken further as shown by Epstein [33] p. 633 who gives, 

3,(4, A- 1 ) = ^ + 7rr 3 C fl (3) + 27rr 3 x(?') + D' X W) , (80) 

where (/' = e -71 "^ and 

X ( g = ^(-l/r)=^ m 3 ■ 

m=l v ' 
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Equation (79) with (80) is appropriate for the high temperature limit since q' 
is the 'inverse' temperature Boltzmann factor. The low temperature form follows 
by simple homogeneous scaling of the Epstein ^-function, which gives, 



Z 2 (2s,A- 1 )=C 2s Z 2 (2s,A), (81) 
and therefore the equivalent form, 

Z 2 (4, A- 1 ) = + ^(3) + 2nC ^(q) + 4^" 2 D X (q) , (82) 

where q = e~^^. 

The low temperature form (82) yields a somewhat simpler expression for the 
free energy, (79), as the Cr(3) term goes out and there is some cancellation that 
produces, 

3 240 2tt AW 

I have thus regained, at some length, the summation form (37) with (39). 

All that has been done is to rederive, in a special case and in a detailed way, 
the standard statistical mode sum, the general form of which was obtained, in this 
fashion, in [34,35] and we are no further forward in finding a closed form for the 
free energy, which is one of my aims. 

The equivalence of (80) and (82), derived here rather trivially (granted Ep- 
stein's calculation), is a known inversion identity, and can also be derived from 
(31). The exact identity is written out in Katayama [36]. His derivation seems 
to be similar to that via the Epstein function. Clearly the higher odd sphere ex- 
pressions will involve the Riemann ^-function at positive odd integers. Katayama 
writes out the one appropriate for the (partial) five-sphere. 

Smart, [16], also derives these identities, which is not surprising since his work 
is concerned with the evaluation of the Epstein ^-function at integer arguments, 
through which he is led to the forms <p 2 k- We see that Epstein had already arrived 
at the same development. Consult also Bodendiek and Halbritter, [37], for a related 
treatment. 
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7. The Epstein-Kober-Selberg-Chowla formula. 



Smart makes use of the, oft quoted, paper of Selberg and Chowla, [38] , which 
is concerned, partly, with the evaluation of the Epstein function via the Kronecker 
limit formula 4 and elegantly yields some of the finite forms for the complete elliptic 
integral used in [2]. The important expression for the Epstein function quoted by 
Selberg and Chowla, [39], and derived by them in [38], (equ.(6)), in terms of K— 
Bessel functions, is essentially the same as that already given by Epstein, [33], 
p. 631, equn.(12), p. 622, equn.(16). Rankin [40] obtained the same result and it also 
occurs in Bateman and Grosswald, [41]. These are the standard references. See also 
Terras, [6] p. 209. Smart, [16], employs this formula to split the Epstein function 
into (simpler) holomorphic and anti-holomorphic parts using the explicit form of 
the Bessel function, K t _ 1 / 2: which henceforth disappears from his analysis. Other 
expressions were given by Deuring [42] and Mordell, [43] . 

An early, and little known, derivation is that by Kober, [44]. Following, and 
generalising, Watson, [45], he obtains, starting from the Bessel function end, the 
formula (I retain his notation), 

i A ( 2w+ i)/ 4 r(^ + i/2) 

01 8tt w+1 / 2 2V ' ' ' ' ' 

u- w T(w) A , n . u w T(w+ 1/2) ^ . 

oo 

o~ 2w (n)n~ w cos(2irvn)K w (27run) , 

n=l 

where the quadratic form in the Epstein function, Z 2 , is amf + 2bm 1 m 2 + cm 2 and 

A i2 _ 2 2 _ ^ 

A = ac — o = a u , v = — . 

a 

Equation (83) is identical to the corresponding one in [33] and [38] and therefore, 
with justice, should be called the Epstein-Kober formula. It is the Fourier expansion 
of the Epstein function in v. For the sum of squares case, v = 0, the result is 
rederived by Guinand [46] in the same way. 

Kober also discusses behaviour under 'reciprocal' transformations, u — > 1/u. 
In particular the diagonal case, v = 0, when one can write u = x. 



4 Landau, [9], gives some interesting history of this famous formula. 
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In this way of doing things, the Bessel function can be eliminated using the 
standard formula 

n s 



E^hr = U s + w )Us-w), (84) 

n=l 

and the Mellin transform (Heaviside's integral), 

'• q _|_ nn q qn POO 

-n-r(-^-)r(-^-) = J dyy'^K^ny), Re a > 1 + |Re W | , (85) 



so that, 

r>c+icso 



E^^^( 2 ™) = t- / ^«- s e(«+^)^-^), (86) 

n=1 n Z7Tl Jc-ioo 



with c > 1 + w and where, 

MEV»rHy2 W ). 

This is the same as Deuring's expression, [42] p. 589. 

The integrand possesses a functional equation which follows most easily from 
that for the Riemann ^-function, £(2w) = £(1 — 2w), and is, 

£(s + u;) ^(a - w) = f w (s) =f w (l-s). (87) 

Using (87) one could derive the functional equation for the Epstein ^-function 
in (83). The usual derivation employs the inversion relation for generalised 9- 
functions, [33,47]. 

Furthermore, it is possible to reverse the argument and use (86) to obtain 
(83) (for v = 0). The diagonal Epstein function is particularly easy to deal with 
being the ^function on the torus, S 1 xS 1 . It can be expressed very simply in terms 
of the (functions on the circle factors (these are Riemann ^-functions ) and the 
result is exactly (83) with (86) without the intervention of the Bessel function. The 
individual Riemann ^-functions in (83) arise from adjusting the zero modes. All 
this is standard and can easily be generalised to higher dimensions. 

A natural step is to proceed as with the holomorphic modular forms and ask 
for the corresponding period polynomials which will follow directly from (86) and 
(87) in the standard manner. 
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There are four poles in the relevant strip, — c < s < c, and a straightforward 
contour calculation produces the result, 

n=l U U n=l U (88) 

= h{2w){u w - 1 - u~ w ) + Ui^w)^-™- 1 - u w ) , 

— ■ — ■ 

The right hand side of (88) might be termed a period function. 

This formula was obtained by Guinand, [46], but using the complete equation 
(83), with v = 0, and properties of the Epstein function. This is unnecessary as I 
have just shown. 

In order to make a connection with the period polynomials encountered ear- 
lier, e.g. (35), one introduces the formula for the K-Bessel function, written as, 

/w^-yf^i-i)'-^—)'"^' * ez - < 89 > 

Then, somewhat similarly to [16], 




so that, for the left-hand side of (88), one finds, 

^MiiT 3 '^ - \ymi~ l St(i/u) ) ■ <9o) 

It is not obvious from (35) that this reduces so nicely to (88), and I leave this as a 
question mark. 

Further expressions involving Bessel functions appear in the Appendix. 
8. Conclusion 

It does not seem possible to ellipticise integrals of the internal energy, in par- 
ticular the free energy (except for the circle), because the modular properties are 
more complicated possessing nonzero cocycle functions. These are related to a 
known Lambert series connected with the Eisenstein series and Zagier's expression 
for the period functions of these is more neatly re-computed by a contour method. 
The use of the fully subtracted series e^ ub allows a more compact treatment. 

Sufficient historical material has been exhibited to indicate that the Selberg- 
Chowla formula should be renamed the Epstein-Kober formula. 
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Appendix, Dirichlet Series. 



It is possible to give a fairly broad context to the inversion behaviour (e.g. 
modular properties), without too much restriction, in terms of Dirichlet series. 
Although these ideas, by now classic, are associated with Hecke, and have been 
expounded by Ogg [48,49], for example, the simplest case occurs earlier in other 
places. I mention Koshliakov, [50], whose discussion I now paraphrase both for 
historical justice and interest. 

Define two ordinary Dirichlet series ('zeta functions'), 

oo oo , 

n n 

n=l n=l 

for some i/'s depending on the asymptotics of a n and b n , and assume that the a n 
and b n are such that the corresponding 'cylinder kernels' satisfy the 'formula of 
transformation' , 

oo oo 

Y. a ^- nbP =^Y. h ^- nb/P i " = "a, P>0, (91) 
n=0 ' n=0 

for fixed a and 6, where the numbers (not necessarily integers) of 'zero modes' are 
taken to be, 

a = -<P(0) b = -iP(0) . 

Then Koshliakov, [50], shows the equivalence of (91) with 

a T(u-s)jj(u-s) _ T(8)4>(8) 

b u ~ s b s ' 

by brute force using the summation function, cr(z), given by the Bessel function 
expression, 

oo , 

a(z) = -2abz^-^ ^ K v _ x (26v^) , 

n=l 

and Heaviside's formula, (85). a(z) has the important property of possessing poles 
at z = 1, 2, . . . with residues a 1; a 2 , . . . and also has a cut along the negative x-axis 
with discontinuity ab u 2ni r ip(Qi)(— x) !/_1 /r(z/), z = x + iy. 

Koshliakov shows that 0(s) is a single-valued holomorphic function having a 
first-order pole at s = v, 

*«~-*< >F£)rb- (93) 
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and it is then a theorem that (92) together with (93) is equivalent to (91), as 
can also be established easily, and instructively, by Mellin transform. It follows, 
symmetrically, that ip(s) is a single-valued holomorphic function having a first- 
order pole at s = v, 

m ~ -0(0) t] , 1 • 

aT(u) s — v 

Koshliakov also gives an integral form of </>(s), 

,. . sin7rs f°° , a(x) n . 
0(s) = / dx^j-, Res<0, (94) 

TV Jq X s 

showing that <j)(s) vanishes at negative integers in agreement with (92). 

Hecke, [15], derives the same equivalence, but only for a n = b n , and eight 
years later. The more general case can also be found in the standard references 
Ogg [48,49] and Weil, [51], and no doubt elsewhere. A related formula involving 
Ramanujan's 'reciprocal functions' occurs in Hardy and Littlewood, [7]. 

Defining cylinder-kernels, or theta-series, including the zero modes, by, 

oo oo 
n=0 n=0 

the reciprocal relation, (91), reads, 

$(6p) = 4*(6/p), (95) 
P 

(continued from the regions of convergence in s) and one has, in standard fashion 
(tf(20)), 

| rOO -j roc 

0(5) = ^^ dp (3 s - 1 -a ), i;(s) = —j^ dp (3 s - 1 (9(0) -bo). 

As an example, the Dirichlet series corresponding to the classical Eisenstein 
series is, from (1) or (84), 

oo 

E a 2t-M) m~ s = ( R (s) ( R (s + 1 - 2t) , 

m=l 

and this is a case discussed by Koshliakov, [50] p. 19, with a n = b n = cr 2t _i(n), v = 
2t, a = (-1)*, b = 2tt, 0(0) = V(0) = ~B 2t /4t. 

Koshliakov also considers Jacobi elliptic cases that yield reciprocal identities 
covered earlier by Glaisher [25]. 
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Generalised Dirichlet series. 



I have described above the classic Dirichlet set up. A certain, and sometimes 
only apparent, generalisation is obtained by replacing n by arbitrary sequences, A m 
and fi n , and defining new and ip by 



OO OO 7 

<k*) = E if > ^) = E # > ( % ) 

m=l 771 n=l ^ n 

which satisfy, by definition, the functional equation (one of many possible), 

T(5 - s)ifj(5 - s) =T(s)<p(s) , 5eR, (97) 
with corresponding heat-kernels, 

oo oo 

m = E a - e-Am/3 ' *(/?) = E 6 » e ~' 1J3 - ( 98 ) 

m=l n=l 

Bochner, [52], (see also Chandrasekharan and Narasimhan, [53], Knopp, [54] ), 
derives, by Mellin transforms, the 'modular relation', 

*(/*)- p-'9(±)=B(P), (99) 
where the 'residual function' B(f3) is 

B(P) = ^~ [ ds X (s) (3~ s = V ^" SRes ^ ■ 

x( s ) is the joint continuation of the two sides of (97) and is assumed to have only 
simple poles for singularities confined to a compact region, 5, of the s-plane, usually 
on the real axis. C(S) is a loop surrounding S. I especially note that if there exist 
higher-order poles then integer powers of log j3 occur. 

In certain situations, e.g. Koshliakov's, B(f3) can be naturally amalgamated 
with the left-hand side to give an exact modular relation like (91), e.g. Bochner, 
[52] p. 342. Typically one would include zero modes A = and fj, = in (98), but 
not in (96). (There seems to be a misprint on p. 342 of [52] which has A = fj, = 1.) 

There is another formula connected with this analysis that is useful. It is 
concerned with the representation of the modified series, 

oo 

4(8, W) = E ( A + W 2). ReW>0 - ( 10Q ) 

n=l ^ n ' 
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There are many reasons why we should wish to add the number w 2 . For 
example it might correspond to a mass or simply be added for extra flexibility as 
when calculating heat-kernel expansion coefficients via resolvents or it might be 
another part of the denominator as in Epstein's derivation of (83); see [55]. 

The representation alluded to is, (Berndt [56]), 

F(s) 0( S , w) = R(s, w) + 2j2bj^) K s _ s (2«^J , (101) 

n=l ^ ' 

where s is such that the summation converges absolutely and R(s,w) is given by 

R(s, w)= T ( s ~ s ') w 2s '~ 2s Res x(s') ■ (102) 
s'es 

Berndt's second proof of (101) proceeds via the cylinder-kernels, $ and I 
prefer his first, Mellin, proof, [57]. Of course, the information employed is the same. 
I give this proof out of interest. 

The Mellin transform relation between 0(s) and the <f>(s,w) of (100), is 

r(s) w) = — ds' w 2s '- 2s T(s - s')T(s') <P(s') 

where c is chosen sufficiently positive to make both <p(s) and ip(s) converge. The 
vertical contour is now moved to the left so as to run from 5 — c — ioo to 5 — c + ioo 
picking up the residues, R(s,w) : (102), on the way. The region, <S, lies between 
these vertical lines and also the horizontal pieces contribute zero. 

On the shifted vertical contour make the coordinate change s' — > 5 — s' and use 
(97) i.e. x($ — s ') = T(s')i(j(s'), valid in this range, and so for this part I find, 

£&n^ / ds' F(s' - 6 + sW^-Jw 25 - 2 *'- 2 * 
= 2 £M^f) K s _ s (2w^Jl n ), 

n=l 



on using (85). Hence (101) has been established. 

In specific situations (101) is quite familiar in the (^function area and in back- 
ground and finite-temperature quantum field theory. 

An illustrative example is the simplest, diagonal Epstein function, 



oo 

v' 1 



(j){s) = Z p (s)= / m m y » m= (m 1 ,m 2 ,...,m p ), (103) 
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which obeys the functional equation (Epstein [33]), 



r( s ) z p ( s ) 



7T 



-^ 2 r(p/2-s)Z p (p/2-s), 



(104) 



so that, see (97), 5 = p/2, a 



1, A 



m.m, b n = tt 



■P/2 



, n n = 7r 2 m.m. 



n 



m 



T(s) Z (s) has poles at s = and s = p/2 and (101) gives for the modified, or 
inhomogeneous, or 'massive' Epstein function, Z Js,w), 



where the dash means to omit the m = term. Note that the w~ 2s could be 
included in the left-hand sum if this were extended to include m = 0. 

The right-hand side can be taken as the continuation of the left-hand side 
showing the single pole at s = p/2, as is correct and may be proved in other ways. 

Actually (105) is shown much more directly using the Jacobi inversion relation, 
which is, of course, how Epstein obtained (104). Equivalent to Jacobi inversion is 
Poisson summation. 

Equations like (105) also occur in lattice summations. Some references, but 
none earlier than those already quoted, are given in the review by Glasser and 
Zucker [58] p. 109. 

A similar formula holds for the general Epstein function which may, or may 
not, have a pole at s = p/2 and may, or may not, vanish at s = 0. The works 
[59,3,60], for example, can be consulted for details, some further references and 
physical applications. 

When p = 1, (f)(s) = 2£ R (2s) and I recover an earlier, well-known formula, 
[45,44], which yields, after putting w = um 2 and summing over m 2 , the formula 
for the Epstein ^-function mentioned earlier; see (83). 

It is left as an exercise to show that the non-diagonal formula, (83), can be 
obtained in the same way using reciprocal relations for the Hurwitz ^-function, or, 
equivalently, the Lipshitz formula ( cf Epstein, [33], Kober [44] p. 622). 

As an application of (105) it might be helpful if I present a standard deriva- 
tion of the usual statistical free energy mode sum from the thermal C^nmction 
expression, for a reasonably general system. The starting point is the determinant 





(105) 
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form, 



oc 



d. 



F = - lim 

2 s^o 



s 



L to i V' 

>/3 s^o s ^ 



(^2 + 47T 2 m 2 //3 2 ) 



n 



s t 



(106) 



where w 2 , and <i n are the eigenvalues and degeneracies of the appropriate operator 
(related to the Laplacian) on the spatial section, Ai. The u> n are the single-particle 
energies. The dash here means that the denominator should never be zero (c/ (74)) 
but the sum includes m = 0. For simplicity it is assumed that there is no spatial 
zero mode and so the dash can be removed. 

Now consider (105) for p = 1. Include the w~ 2s with the left-hand sum, set 
w = (3uJ n /27r, multiply by the degeneracy, d n and sum over n. The limit (106) is 
then easily taken since sT(s) = T(s + 1) and only the K 1 / 2 Bessel appears. There 
is nothing deep about this calculation and the result is a standard form, 



To avoid specifically field theoretic problems, it has been assumed that ( M (— 1/2) 
exists. This is true for the cases discussed in this paper. 

The result, (107), is an analogue of the (first) Kronecker limit-formula applied 
to the generalised torus S 1 x M; see [61]. The original formula relates to M =S X 
and, in this case, there exists a functional relation, (97), so that the limit-formula 
is often expressed in terms of the remainder about the pole in ((s, j3) at s = 1. In 
the general case this is not possible. 

Finally, if one introduces the standard arithmetic quantity r p (n), equal to the 
number of representations of the integer n as the sum of p squares, then the diagonal 
Epstein function, (103), is obtained by writing, 



and is treated as such by Koshliakov who also gives the 'exact', product forms of 




(107) 



^C^(-l/2) + ^rf n log(l-e-^). 



n 




n 



Z p (s) for p = 2 and p = 4. 
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